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Abstract 

Q^' We study the finite-temperature effective potential of the Minimal Supersymmetric 

Standard Model in the full (mA,tan/9) parameter space. As for the features of the 
electroweak phase transition, we identify two possible sources of significant differ- 
ences with respect to the Standard Model: a stop sector with little supersymmetry 
breaking makes the phase transition more strongly first-order, whereas a light CP- 
5^ I odd neutral boson weakens its first-order nature. After including the leading plasma 

effects, T = radiative corrections due to top and stop loops, and the most impor- 
tant experimental constraints, we find that the danger of washing out any baryon 
asymmetry created at the electroweak scale is in general no less than in the Standard 
Model. 
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1. The option of generating the cosmological baryon asymmetry at the electroweak 
phase transition is not necessarily the one chosen by Nature, but it is certainly fascinat- 
ing, and has recently deserved a lot of attention [|l|. At the qualitative level, the Standard 
Model (SM) meets the basic requirements for a successful implementation of this scenario. 
At the quantitative level, however, it suffers from two basic problems. First, the amount 
of CP violation in the Kobayashi-Maskawa matrix appears to be insufficient, even after 
taking into account the large uncertainties associated with dynamical models of baryogen- 
esi^. Secondly, one must make sure that sphaleron interactions in the broken phase do not 
wash out, at the completion of the phase transition, any previously created baryon asym- 
metry. This requirement suggests that the transition should be rather strongly first order, 
v{Tc)/Tc ^ 1, where Tc is the critical temperature and v{Tc) is the symmetry-breaking 
vacuum expectation value. In the SM, this condition turns out to be incompatible with 
the experimental limits 0] on the Higgs mass, > 63.5 GeV at 95% c.l., even after 
implementing the conventional techniques for dealing with the infrared problem To 
cure both problems, one can consider plausible extensions of the SM, to see if they can 
allow for additional sources of CP violation and for an enhanced strength of the first-order 
phase transition. Among these extensions, the physically most motivated and phenomeno- 
logically most acceptable one is the Minimal Supersymmetric Standard Model (MSSM). 
This model allows for extra CP-violating phases besides the Kobayashi-Maskawa one, 
which could help in generating the observed baryon asymmetry 0. It is then interesting 
to study whether in the MSSM the nature of the phase transition can be significantly 
modified with respect to the SM. 

In a recent paper [0, some of us have considered the MSSM in the limit itla oo, 
corresponding to only one light Higgs with SM properties, and improved over previous 
studies PI by including a full discussion of the top/stop sector, and by resumming the 
leading plasma corrections to gauge boson and stop masses. It was found that this special 
limit of the MSSM can only marginally improve the situation with respect to the SM 
case. In the present paper, we extend the considerations of to the full (m^,tan/5) 
parameter space, characterizing the Higgs sector of the MSSM. Even barring the interesting 
possibility of spontaneous CP- violation at finite temperature as well as the possibility 
of charge- and colour-breaking minima, we have to deal with a complicated two-variable 
potential, which requires a numerical analysis. However, to allow for an understanding of 
the behaviour of various quantities, we also produce some approximate analytical formulae. 
After including the most important experimental constraints, we find that there is very 
little room for the MSSM to improve over the SM. 

2. The main tool for our study is the one-loop, daisy-improved finite-temperature 
effective potential of the MSSM, Vcs{4>,T). We are actually interested in the dependence 
of the potential on (pi = KeH^ and 02 = Reifg o'^ly? where and if^ the neutral 

^Unorthodox views on this point, recently put forward in have been subsequently questioned in 
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components of the Higgs doublets Hi and H2, thus (f) will stand for 
the 't Hooft-Landau gauge and in the DR-scheme, we can write 



h,4>2)- Working in 



Kfr(0, T) = Vo{(P) + Vi{(P, 0) + AVi{(P, T) + AKiaisy(0, T) 
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The four contributions (IH^) to the effective potential (|l|) have the following meaning. The 
first term, eq. (||), is the tree-level potential. The second term, eq. (|]), is the one-loop 
contribution at T = 0: Q is the renormalization scale, where we choose for definiteness 
Q"^ = m|, mf{(t)) is the field-dependent mass of the i^^ particle, and rii is the corresponding 
number of degrees of freedom, taken negative for fermions. Since Vi((/>, 0) is dominated by 
top it) and stop (^1,^2) contributions, only these will be included in the following. The 
third term, eq. (|), is the additional one- loop contribution due to temperature effects. 
Here Jj = J+(J_) if the i^^ particle is a boson (fermion), and 



dx X log 1 =F e 
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Since the relevant contributions to AVi(0, T) are due to top (t), stops (ti,t2) and gauge 
bosons {W, Z), only these will be considered in the following. Finally, the last term, eq. (^), 
is a correction coming from the resummation of the leading infrared-dominated higher-loop 
contributions, associated with the so-called daisy diagrams. The sum runs over bosons 
only. The masses m^(0, T) are obtained from the m'f{(f)) by adding the leading T-dependent 
self-energy contributions, which are proportional to T^. We recall that, in the gauge boson 
sector, only the longitudinal components {Wl, Z^^'-^i) receive such contributions. 
The relevant degrees of freedom for our calculation are: 
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The field-dependent top mass is 



The entries of the field-dependent stop mass matrix are 



ml^+ml{ct^) + Dl{cl^) 
,2 I _2/j,\ I n2 



(9) 
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where mq^, mu^ and At are soft supersymmetry-breaking mass parameters, is a super- 
potential Higgs mass term, and 
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are the D-term contributions. The field-dependent stop masses are then 



ml {4>) + m}U) 
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(14) 



The corresponding effective T-dependent masses, mf^ ^{(j),T), are given by expressions 
identical to (|14|) , apart from the replacement 

(15) 



The Ilr (T) are the leading parts of the T-dependent self-energies of n , 
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where gs is the strong gauge coupling constant. Only loops of gauge bosons, Higgs bosons 
and third generation squarks have been included, implicitly assuming that all remaining 
supersymmetric particles are heavy and decouple. If some of these are also light, the 
plasma masses for the stops will be even larger, further suppressing the effects of the asso- 
ciated cubic terms, and therefore weakening the first-order nature of the phase transition. 
Finally, the field-dependent gauge boson masses are 



2' ' ^' ' 2 

and the effective T-dependent masses of the longitudinal gauge bosons are 
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In eqs. (|r^) and (^), IIwl(T) and IIbj^(T) are the leading parts of the T-dependent 
self-energies of Wl and B^, given by 



^H^.(T) = ^/T^ ^^JT) = ^^'^T^ 



1^ 
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where only loops of Higgs bosons, gauge bosons, Standard Model fermions and third- 
generation squarks have been included. 



3. We shall now analyse the effective potential (|I|) as a function of and T. Be- 
fore doing this, however, we trade the parameters 1711,7712, ml appearing in the tree-level 
potential (|^) for more convenient parameters. To this purpose, we first minimize the zero- 
temperature effective potential, i.e. we impose the vanishing of the first derivatives of 
Vo(0) + Vi(0, 0) at (0i,</>2) = {vi,V2), where (^1,^2) are the one-loop vacuum expectation 
values at T = 0. This allows us to eliminate mf and in favour of m| and tan f3 = V2/V1: 
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Moreover, can be traded for the one-loop-corrected mass m\ of the CP-odd neutral 
Higgs boson. In our approximation 
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Therefore the whole effective potential (||) is completely determined, in our approximation, 
by the parameters (m^,tan/3) of the Higgs sector, and by the parameters (m^, mq^, m^/g, 
/i. At) of the top/stop sector. The same set of parameters also determines the one-loop- 
corrected masses and couplings of the MSSM Higgs bosons. 

The next steps are the computation of the critical temperature and of the location of 
the minimum of the effective potential at the critical temperature. We define here Tq as 
the temperature at which the determinant of the second derivatives of Ves{4',T) at </> = 
vanishes: 
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It is straightforward to compute the derivatives in eq. (^Sf) from the previous formulae; 
the explicit expressions are 
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The coefficients aij,bij are given by 
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and the coefficients aij,bij are given by identical expressions, apart from the replacement 
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Once eq. (|25|) is solved (numerically) and Tq is found, one can minimize (numerically) 
the potential V^siij^.T^ and find the minimum [t'i(To), f2(To)]. The quantity of interest is 
indeed, as will be discussed later, the ratio f (To)/To, where f (Tq) = \Jv'1{Tq) + V2{Tq). 



4. Before moving to the discussion of our numerical results, we would like to present 
some approximate analytical formulae, which will be useful for a qualitative understanding 
of the various dependences of Tq, fi(To) and f2(To). In this paragraph, we shall work in 
the limit of heavy degenerate stops, mgg = mjj^ = m ^ Tq and = /i = 0, neglecting 
the D-term contributions to the stop squark masses, and keeping only the most important 
terms in the high-temperature expansions of the Jj functions for the gauge bosons and 
the top quark. In the chosen limit, the effective potential of eq. (|1]) can be approximately 
written, in the polar coordinates = y 0i + 01 tan 6' = 02/0i, as 



Kff(0, T) ~ [a{e)T^ - h{e)\ - E V ^ 1 At(^)^' , (29) 
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We have exploited the fact that, in the high-temperature limit, the terms proportional to 
mf log cancel in the sum Vi{(j), 0) + AVi{(j), T) (in the chosen limit, of course, we cannot 
perform the high-temperature expansion on the stop contributions). 

We define the critical angle 6* by the flat direction of the effective potential (|I]) around 
= at T = To, and we denote by 9^^^ the analogous quantity evaluated from the 
approximate parametrization of eq. (^91) . The critical temperature Tq and the critical 
angle 6'*pp are then determined by the conditions 
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Solving eq. (|34) amounts to solving eq. (^Sf) and finding the eigenvector corresponding to 
the zero eigenvalue. One finds 
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where Tq is determined by eq. (|35|). If we now assume that tan0(To) = V2{Tq) / vi{To) can 
be approximated by tan6'*pp, we can also write 



viTo) 



app -^ro(^app 



. (40) 



In fig. 1, we plot tan^(To) (solid lines) and tan^* (dashed lines) as functions of itla, 
for tan/3 = 1,2,5 and the representative parameter choice rrit = 150 GeV, ttlq.^ = mu^ = 
1 TeV, At = fi = 0. We can see that, in all cases, tan6'(To) ~ tan6'* to a very good 
accuracy. To check our analytical approximation, we also plot tan6'*pp (dotted lines) 
as obtained from eq. (pQ)). We can see from fig. 1 that all three quantities tend to the 
corresponding value of tan/3 for large values of m^, whereas they increase to large values 
for small values of m^. This fact is a general trend of V^fj, for arbitrary values of the 
top/stop parameters. As for f (To)/To, we have also checked that the analytical expression 
(|40| ) is an adequate approximation to the numerical value obtained from (|l]) and (^). 
The qualitative behaviour of v(Tq)/Tq can then be derived from (0), (^), (^) and (^): 
for fixed tua, v{Tq)/Tq increases when tan/3 is approaching 1 from above; for fixed tan/3, 
v{Tq)/Tq is an increasing function of uia- 



5. We now discuss the particle physics constraints on the parameters of the top/stop 
sector and of the Higgs sector. To be as general as possible, we treat m^g, m^/g and the 
other soft mass terms as independent parameters, even if they can be related in specific 
supergravity models. 

The constraints on the top/stop sector have already been discussed in [0, so we just 
recall them briefiy. Direct and indirect searches at LEP 0] imply that <; 45 GeV, 
which in turn translates into a bound in the (mgg, tan /3) plane. Electroweak precision 
measurements [|rT| put stringent constraints on a light stop-sbottom sector: in first ap- 
proximation, and taking into account possible effects |T2| of other light particles of the 
MSSM, we conservatively summarize the constraints by Ap{t,b) + Ap{t,b) < 0.01, where 
the explicit expression for Ap(t,b) can be found in [p!3| . 

We finally need to consider the constraints coming from LEP searches for supersym- 
metric Higgs bosons [Q. Experimentalists put limits on the processes Z hZ* and 
Z — > hA, where h is the lighter neutral CP-even boson. We need to translate these 
limits into exclusion contours in the (mA,tan/3) plane, for given values of the top/stop 
parameters. In order to do this, we identify the value of BR{Z hZ*), which corre- 
sponds to the limit > 63.5 GeV on the SM Higgs, and the value of BR{Z —>■ hA), 
which best fits the published limits for the representative parameter choice rrit = 140 GeV, 
= mif^ = m = 1 TeV, At = fi = 0. We then compare those values of BR{Z hZ*) 
and BR{Z hA) with the theoretical predictions of the MSSM, for any desired parameter 
choice and after including the radiative corrections associated to top/stop loops [Il4| , |l0 



Of course, this procedure is not entirely correct, since it ignores the variations of the ef- 
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ficiencies with the Higgs masses and branching ratios, as well as the possible presence of 
candidate events at some mass values, but it is adequate for our purposes. 



6. We now present our numerical results, based on the effective potential of eq. (|l|), 
concerning the strength of the electroweak phase transition and the condition for preserving 
the baryon asymmetry. According to [|15|, the condition to avoid erasing any previously 
generated baryon asymmetry via sphaleron transitions is 

Es^H^ > 45 , (41) 
Tc 

where Tq is the actual temperature at which the phase transition occurs, satisfying the 
inequalities 

To<Tc<Tn, (42) 

if To is defined by ( |25| ) and Tjj is the temperature at which there are two degenerate 
minima. Particularizing to the MSSM the studies of sphalerons in general two-Higgs 
models |T^, we obtain that 

<r'^(T)<i?S)f(T), (43) 

where, in our conventions, 

gg'lr) 4V2^ f \^le(T)] \ v{T) 



and 5 is a smoothly varying function whose values can be found in |T^. For example 



5(10-2) = 1.67, 5(10-1) = 1.83, B{1) = 2.10. It can also be shown that 

vjTnl^viTcl^viT^ 
Td Tc To ' ^ ^ 

Finally, the corrections in due to g' have been estimated and shown to be small 
|1T8| . Therefore, a conservative bound to be imposed is 

R^-^ ^ >1- (46) 

To 4:5g 

The last point to be discussed is the determination of the value of AefT[^^(To)] to be 
plugged into eq. (^). The i?-function we use is taken from ref. [T^, where the sphaleron 



energy was computed using the zero-temperature 'Mexican-hat' potential, V = ^{(p'^—v'^^'^ 



The sphaleron energy at finite temperature was computed in ref. [|TU| , where it was proven 
that it scales like v(T), i.e. 



i^3f(T)=i^J^(0)^, (47) 
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with great accuracy. Therefore, to determine the value of Aefr[^(To)] we have fitted 
Vcs{4>,To), as given by eq. (P, to the appropriate approximate expression, 

Ves{(p, To) ~ iAcfr[^(To)][0^ - f ^(Tq)]^ + field-independent terms , (48) 
where the field-independent terms are just to take care of the different normalizations of 



the left- and right-hand sides. The value of Aefr obtained from (^Sj), 



where all quantities on the right-hand side are calculated numerically from the potential 
of eq. dll), is then plugged into eq. (|^) to obtain our bounds. We have explicitly checked 



the quality of the fit in eq. (^81) , finding an agreement that is more than adequate for our 
purposes. 

Our numerical results are summarized in fig. 2, in the (mA,tan/3) plane and for two 
representative values of the top quark mass: rrit = 130 GeV (fig. 2a) and = 170 GeV 
(fig. 2b). In each case, the values of the remaining free parameters have been chosen in 
order to maximize the strength of the phase transition, given the experimental constraints 
on the top-stop sector. Notice that arbitrarily small values of mc/g cannot be excluded 
on general grounds, even if they are disfavoured by model calculations. Also, we have 
explicitly checked that, as in ref. 0, mixing effects in the stop mass matrix always worsen 
the case. In fig. 2, solid lines correspond to contours of constant R: one can see that the 
requirement of large values of R favours small tan/3 and ^ m^. The thick solid line 
corresponds to the limits coming from Higgs searches at LEP: for our parameter choices, 
the allowed regions correspond to large tan (3 and/or itla ^ mz- For reference, contours 
of constant rrih (in GeV) have also been plotted as dashed lines. One can see that, even 
for third-generation squarks as light as allowed by all phenomenological constraints, only 
a very small globally allowed region can exist in the (m^, tan /3) plane, and that the most 
favourable situation is the one already discussed in ref. More precisely, the region 
that is still marginally allowed corresponds to rriA S> mz, tan/3 ~ 2, stop and sbottom 
sectors as light as otherwise allowed, a heavy top, and a light Higgs boson with SM- 
like properties and mass rrih ~ 65 GeV, just above the present experimental limit. A less 
conservative interpretation of the limits from precision measurements, the inclusion of some 
theoretically motivated constraints on the model parameters, or a few GeV improvement in 
the SM Higgs mass limit, would each be enough to fully exclude electroweak baryogenesis 
in the MSSM. 



7. In summary, our analysis of the full (m^, tan /?) parameter space extends and 
confirms the results of ref. in the region of the MSSM parameter space allowed by 
the present experimental constraints, there is very little room for fulfilling the constraint 
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(|46|), which is a necessary condition for electroweak baryogenesis. To put these results in 
a clearer perspective, some final comments on possible ways out are in order. 

First, one could think of relaxing the constraint tan/3 > 1 (and the corresponding LEP 
bounds), which is usually motivated by the theoretical assumption of universal soft Higgs 
masses at the SUSY-GUT scale, Mu ~ 10^^ GeV. The possibility of tan/? < 1, however, is 
incompatible with a heavy top quark, since, for > 130 GeV and supersymmetric particle 
masses of order mz, the running top Yukawa coupling would become non-perturbative 
at scales smaller than Mu: such a possibility is strongly disfavoured by the successful 
predictions of the low-energy gauge couplings in SUSY GUTs. 

A second possibility is that large non-perturbative effects, neglected by conventional 
calculational techniques, modify the predicted values of the sphaleron energy and/or of 
v{To)/To (for recent suggestions along this line, see [0). We do not see strong physical 
arguments to favour this, but we admit that it cannot be rigorously excluded. Perhaps 
alternative approaches to the electroweak phase transition could help clarify this point 
in the future. 

Barring the above-mentioned possibilities, one could still try to rescue electroweak 
baryogenesis by further enlarging the MSSM Higgs sector, for example by introducing an 
extra singlet. Supersymmetric models with singlets and non-supersymmetric models, how- 
ever, develop dangerous instabilities if coupled to the superheavy sector of an underlying 
unified theory. It might well be that baryogenesis has to be described by physics at a scale 
larger than the electroweak one. 
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Figure captions 



Fig.l: The quantities tan 9{To) (solid lines), tan 9* (dashed lines) and tan 6'*pp (dotted lines), 
as functions of m^, for tan /? = 1, 2, 5 and the representative parameter choice rrit — 
150 GeV, = mu^ ^ I TeV, At = n^O. 

Fig. 2: Contours of R in the (m^, tan/3) plane, for the parameter choices: a) rrit = 130 GeV, 
mq, = 50 GeV, mu, = (m^- ~ 130 GeV, m^^ ~ 50 GeV), At = fi = 0; 
b) mt = 170 GeV, tuq.^ = 280 GeV, = (m^^ ~ 330 GeV, rrii^ ~ 170 GeV, 
m^^ ~ 280 GeV), At = /J, — 0. The region excluded by Higgs searches at LEP is 
delimited by the thick solid line. For reference, contours of constant ruh (in GeV) 
are also represented as dashed lines. 
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